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We present results, for heavy to moderate quark masses, of a study of thermodynamic properties of 1-flavour 
QCD, using the multiboson algorithm. Finite-size scaling behaviour is studied on lattices of size 8 3 x 4, 12 3 x 4 
and 16 3 x 4. It is shown that, for heavy quarks, the peak of the Polyakov loop susceptibility grows linearly with 
the spatial volume, indicating a first order phase transition. The deconfinement ratio and the distribution of the 
norm of the Polyakov loop corroborate this result. For moderately heavy quarks the first-order transition weakens 
and becomes a crossover. We estimate the end point of the first-order phase transition to occur at a quark mass 
of about 1.6 GeV. 



1. Introduction 

We use the non-hermitian variant of the multi- 
boson method [jlj to simulate one-flavour QCD 
at finite temperature (for more details on the al- 
gorithm see [|[[§). We- perform a finite-size scal- 
ing study for four k values of Wilson fermions. 
Because we simulate heavy quarks only, we can 
safely ignore the problem of the sign of the 
fermionic determinant. Global observables, such 
as the Polyakov loop, obtained with our multi- 
boson algorithm showed an autocorrelation time 
O(10) shorter than if obtained with an efficient 
polynomial variant of HMC, appropriate to sim- 
ulate 1 flavour Q . This substantial improvement 
can presumably be attributed to the additional 
freedom due to the bosonic fields of the multibo- 
son algorithm enabling movement around energy 
barriers in constrast to HMC where one has to go 
over the barrier. 

In Table 1 we collect the parameters of our sim- 
ulations. 
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Table 1 

Summary of simulations: for each n value and vol- 
ume studied, Nb indicates the number of bosonic 
fields used and acc the average acceptance; Ksw 
is the total number of thermalized configurations 
(in thousands) used in the reweighting procedure. 





8 3 x 4 12 3 x 4 16 3 x 4 




Nb/acc Ksw 


Nb/acc Ksw 


N},/ acc Ksw 


0.05 


8/0.78 18 


12/0.74 20 


24/0.83 20 


0.10 


16/0.67 45 


24/0.63 50 


32/0.67 27* 


0.12 


24/0.74 55 


32/0.67 40 


40/0.69 10* 


0.14 


32/0.77 60 


40/0.70 37 


50/0.67 12 



* indicates that we are still accumulating statistics. 



2. Determination of (3 C 

To identify the critical line (3(k) we study the 
signals for a first order phase transition on a finite 
lattice (e.g. ||) using three observables: 

• The histogram of the norm of the Polyakov 
loop, which shows 2 peaks of equal area 



/c=0.05 



at criticality. For the lattices under study 
here we observed enough tunneling events 
between the two phases to make this a re- 
liable method. This is shown in Fig. 1 for 
K = 0.05 for the three lattice sizes. 

• The deconfmement ratio p = 3/2p — 1/2, 
where p is the probability for the complex 
trace of the Polyakov loop to fall within 
±20deg of a Z3 axis, p goes from (uni- 
form angular distribution) to 1 (alignment 
with Z3 axis). 

• The peak value of the susceptibility of the 
Polyakov loop. 

Using these criteria we determined the critical 
values (3(k) given in Table 2. 



K 


0.05 


0.10 


0.12 


0.14 


Pc 


5.692 


5.66 


5.63 


5.59 



Table 2 

The critical values (3{k). The estimated error is 
about one on the last digit. 

3. Order of transition 

In Fig. 2 we show the deconfmement ratio p 
obtained using reweighting q] for k = 0.05, k = 
0.10 and n = 0.12. Across the transition region 
the slope of p((3) increases with the volume for 
k = 0.05 signaling a first order transition. This 
is to be contrasted with the behaviour of p for 
k = 0.12 indicating that k = 0.12 is already in 
the crossover region. For k = 0.10 we observe 
a qualitative behaviour more appropriate for a 
weak first order transition. 

In Fig. 3 we display the peak of the susceptibil- 
ity versus the spatial volume V. The lines shown 
are best fits to the form V a . For n = 0.05 the 
best fit yields a = 0.96(4) whereas for k = 0.14 
a is zero. For k = 0.10 and n = 0.12 the situa- 
tion is less clear. The small value of a = 0.24(4) 
at k = 0.12 as well as the absence of tunneling 
indicate that we are in the crossover region. For 
k = 0.10 a = 0.54(3), increasing with the statis- 
tics of the L — 16 lattice (still in progress), and 
tunneling is still observed. If this increasing trend 
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Figure 1. Tunneling between the two phases for 
k = 0.05. By |f2| we denote the norm of the 
spatial average of the Polyakov line. 

of the peak of the susceptibility for L = 16 con- 
tinues then we will have confirmed that the tran- 
sition is still first order. 

Considering the signals that we have for k — 
0.1, namely tunneling between the two phases, 
the behaviour of the deconfmement ratio and the 
finite size scaling of the susceptibility, we are led 
to the conclusion that the end point of the first 
order phase transition occurs very near n = 0.1. 

Taking the end-point value of k, n ep as 0.1 ± 
0.01, we can approximately map to physical 
units, using the tadpole-improvement property 
k c (^)(D) 1 / 4 « 1/8 to obtain m q a - 1.8(5), and 
(4a) _1 ~ 220MeV from the deconfmement tem- 
perature. This gives m q ~ 1.6(5)GeV at the end- 
point. 
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Figure 2. Deconfincment ratio for k = 0.05, k = 
0.10 and n = 0.12 for three lattice sizes 

This is in line with phenomenological expecta- 
tions. The pure gauge deconfincment transition 
is fairly weak, with a correlation length 0{& few 
0.-1/2^ This is the minimum system size nec- 
essary to observe the deconfinement transition. 
Dynamical quarks introduce a new length scale, 
the distance where the string breaks, 0(2m q /a). 
Confinement can only be observed up to this dis- 
tance. When the quark mass is lowered suffi- 
ciently that the second length-scale is similar to 
(or smaller than) the first, one cannot tell if the 
system is confined or deconfined, and the transi- 
tion is replaced by a crossover. This occurs for 
Tfiq ~ 0(afew^/a/2), i.e. m q ~ 0(l)GeV. 
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Figure 3. Volume dependence of the maximum of 
the Polyakov loop susceptibility 
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